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Height  variations  along  sea  ice  pressure  ridges  were  studied  to  determine 
whether  low  spots  in  the  ridges  persisted  over  long  enough  distances  for  a 
wide  vehicle  such  as  the  surface  effect  vehicle  to  pass  through  them.  Heights 
along  eight  pressure  ridges  varying  in  length  from  0.9  km  to  2.2  km  were 
measured  at  5-m  intervals  using  aerial  photographic  data  obtained  over  the 
Beaufort  Sea.  The  probabilities  of  finding  holes  at  the  mean  heights  of  the 
r  dges  and  0.67  m  above  and  below  the  mean  heights  were  calculated  as  a  func¬ 
tion  of  the  hole  widths  using  an  autocorrelation  function  obtained  from  the 
leight  data.  The  curves  obtained  were  in  good  agreement  with  passage 
probabilities  calculated  directly  from  the  data.  The  results  indicate  that 
the  probability  of  finding  passage  through  a  hole  at  least  0.67  m  below  the 
mean  ridge  height  is  less  than  0.1  for  holes  greater  than  10  ra  wide.  It  is 
concluded  that  a  large  vehicle  such  as  a  surface  effect  vehicle  teat  cannot 
cross  ridges  at  their  mean  heights  will  have  to  go  around  such  ridges  because 
low,  wide  holes  do  not  appear  with  sufficiently  high  probability  to  make  hole- 
-  'arching  a  useful  method. 


14.  Key  Words 

Aerial  photography 
Ice  navigation 


M0 


Pressure  ridges 
Sea  ice 


FN>  4  A  O  •■Pkicia  oo  >mh  i<ii  i  jam  i«.  mien 
IW«  III#  /O  OIMLITI  FOD  llwr  ull 


"I  AT  If  “J  A‘ 
IN f  Oh’  .'A  I 


t.i  irjifAt 
.  '  FRVV 


Surface  effect  vehicles 


Unclassified 

Security  CUiii/lcttion 


HEIGHT  VARIATION  ALONG  SEA  ICE 
PRESSURE  RIDGES  AND  THE  PROBABILITY  OF 
FINDING  “HOLES"  FOR  VEHICLE  CROSSINGS 


W.D.  Hibler  III  and  S.F.  Ackley 


December  1973 


PREPARED  for 

ADVANCED  RESEARCH  PROJECTS  AGENCY 
ARP  A  UkDER  1615 

BY 

(  °RPS  OF  ENGINEERS,  U.S.  ARMY 

COLD  REGIONS  KESEARCH  AND  ENGINEERING  LABORATORY 

HANOVER,  NEW  HAMPSHIRE 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


II 


PREFACE 

This  report  was  prepared  l>v  Dr.  W.D,  Hibler  III.  Research  Physicist,  and  Mr.  S.E. 
Ackley,  Research  Physicist,  of  the  Snow  and  Ice  Branch,  Research  Division,  U.S.  Army 
Cold  Regions  Research  and  Engineering  Laboratory  (USA  CRREL). 

The  study  was  conducted  in  support  of  the  Advanced  Research  Projects  Agency 
(ARPA)  Arctic  Surface  Effect  Vehicle  Program  under  ARPA  Order  1615. 

This  report  was  technically  reviewed  by  Dr.  R.A.  Liston  and  Dr.  Y.  Nakano  of 
USA  CRREL. 

The  authors  wish  to  thank  Mr.  S.J.  Mock  'or  assistance  in  computer  programming  and 
Messrs.  A.  Hartwell  and  J.  Riordar  for  assistance  in  reduction  of  the  ridge  profile  data. 
They  also  wish  to  acknowledge  the  aircraft  support  of  the  NASA  CV-990  (Galileo)  in  this 
study. 

The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  pro¬ 
motional  purposes.  Citation  of  trade  nanxcs  does  not  constitute  an  official  endorsement 
or  approval  of  the  use  of  such  commercial  products. 


Manuscript  received  14  September  1973 


CONTENTS 


Page 

Introduction .  1 

Experimental  data .  3 

Theory .  3 

Results .  6 

Conclusions .  8 

Literature  cited .  9 

Abstract  . 11 

ILLUSTRATIONS 

Figure 

1.  Oblique  aerial  view  of  pressure  ridges  m  the  Beaufort  Sea .  1 

2.  Aerial  photo  triangle  taken  at  1007-in  altitude .  2 

3.  Digitized  profiles  at  friii  intervals  of  ridges  m  Figure  2 .  1 

•1.  Average  probability  of  finding  passage  versus  passage  width .  7 

5.  Percentage  of  fiO-m  ridge  segments  containing  holes  of  various  widths. ...  8 

TABLES 

I.  Mean  height s  and  variances  for  ridges  2-9 .  7 

II.  Average  autocorrelation  function  C(r)  for  ridges  2-9 .  7 


HEIGHT  VARIATION  ALONG  SEA  ICE  PRESSURE  RIDGES  AND  THE 
PROBABILITY  OF  FINDING  “HOLES”  FOR  VEHICLE  CROSSINGS 


by 

W.D.  Milder  Hi  and  S.I\  Ackley 


Introduction 

A  significant  parameter  required  for  the  design  of  a  surface  vehicle  to  traverse  the  sea  ice  of 
the  Arctic  Basin  is  an  estimate  of  its  ability  to  traverse  the  nn  jor  terrain  obstacles,  sea  ice 
pressure  ridges.  This  problem  ts  postulated:  If  a  portion  of  a  ridge  with  points  or  segments  higher 
than  the  vehicle’s  traversing  capability  is  encountered,  bow  far  along  that  ridge  would  the  vehicle 
have  to  travel  to  encounter  a  hole  wide  enough  and  low  enough  to  go  through?  One  possible  solu¬ 
tion  ts  to  travel  until  tin*  ridge  dies  off.  Figure  1,  which  gives  an  oblique  aerial  view  of  first-year 
pressure  ridges  m  a  segment  of  the  Beaufort  Sea,  shows  that  this  is  a  practical  solution  only  it  it 
is  possible  to  cross  most  ridges  either  by  going  through  them  or  over  them.  Clearly,  because  of 
the  networks  formed  by  tin*  ridges,  it  is  nearly  impossible  to  travel  between  two  points  without 
crossing  some  ridges.  Given  the  need  for  the  vehicle  to  cross  ridges,  an  estimate  of  tin*  expecta¬ 
tion  of  holes  of  various  heights  and  widths  existing  in  the  ridges  is  then  required  for  optimum 
vehicle  design.  For  this  reason,  the  height  variations  along  sea  i  *e  pressure  ridges  were  studied 
to  determine  whether  low  spots  persisted  over  long  enough  distances  for  a  wide  vehicle  such  as  the 
surface  effect  vehicle  to  pass  through. 


Figure  1.  Oblique  aerial  view  of  pressure  ridges  in  the  Beaufort  Sea. 


Figure  2.  Aerial  photo  triangle  token  at  1067-m  altitude.  a.  ft,  and  BC  refer  to  strum  measurement  line  end  points.  Numbered  boxes  shou  the  pres¬ 
sure  ridge  locations  investigated  for  height  variations.  Ground  truth  was  obtained  on  ridgt  s  1  and  S  in  this  figure. 
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Experimental  Data 

To  obtain  height  variations  over  sufficient  ridge  lengths  for  use  in  the  statistical  calculation 
shown  in  the  next  section,  aerial  photography  was  used.  This  method  was  used  liecause  ground- 
based  survey  measurements  were  insufficient  and  obtaining  them  would  be  too  time  consuming  m 
a  future  field  program.  However,  using  aerial  photography  to  obtain  the  height  variations  also 
preseinad  a  problem.  To  resolve  small  height  variations  (about  0.2  m),  the  aerial  cameras  must  be 
flown  at  low  level  (500-m  to  1200-m  altitudes).  In  the.  existing  photos  taken  at  low  level,  the 
overlap  between  adjacent  photos  is  usually  only  10  to  15%,  which  does  not  provide  adequate 
stereo  coverage  to  measure  vertical  relief  by  photogrammetne  methods. 

To  overcome  this  problem,  the  technique  of  measuring  shadow  profiles  of  ridges  to  obtain  ridge 
height  measurements  was  employed.  Figure  2  shows  three  strips  of  sea  ice  photos  of  the  strain 
triangle  legs  (Hibler  et  al.  1973)  measured  during  the  AIDJEX  Program  in  1971.  These  photos  were 
obtained  by  RC-8  aerial  camera  flown  by  the  NASA  CV-990  at  an  altitude  or  10(37  m  (3500  ft)  on 
11  March  1971.  The  pressure  ridges  can  clearly  be  seen  in  these,  photos  and  upon  close  inspection 
the  shadows  of  the  ridges  also  can  clearly  he  seen.  Tin*  sun  angle  at  this  time  of  year  is  quite  low, 
so  the  ridges  cast  a  long  shadow.  From  the  altitude  of  the  photography,  the  date,  time  of  day  and 
latitude  of  the  area,  these  horizontal  shadow  lengths  can  be  calibrated  to  give  the  vertical  height 
of  adjacent  points  on  the  ridge.  The  shadow  length  £,  from  a  vertical  object  of  height  h  is 


tan  « 

where  «  is  the  sun  angle  above  the  horizon.  By  measuring  t lie  variations  in  shadow  lengths  along 
a  given  ridge,  the  vertical  height  distribution  can  therefore  lie  computed.  This  method  does  not 
require  stereo  aerial  photo  coverage  since  the  measurement  is  in  the  horizontal  plane  (the  shadow 
length)  and  is  not  directly  of  the  vertical  height  of  the  ridge. 

The  original  9-in.  x  9-in.  photos  were  enlarged  to  approximately  three  times  their  original  size 
to  obtain  a  convenient  working  scale.  Comparisons  of  two  profiles  obtained  by  measuring  the 
shadow  profiles  of  ridges  with  two  profiles  obtained  from  ground  surveys  indicated  an  absolute  error 
o(  approximately  *  0.25  m  in  any  given  height .  The.  relative  error  in  heights  between  any  two  adja¬ 
cent  shadows  was  less  than  this,  of  the  order  of  '  0.10  m.  A  shadow  length  was  obtained  every  5  m 
along  the  ridges  indicated  in  Figure  2.  The  ridges  were  "first-year”  ice  ridges,  with  the  exception 
of  ridge  1  which  was  a  multiyear  ridge.  Details  of  the  structure  of  this  ridge  me  presented  by 
Kovacs  et  al.  (1972). 

The  height  profiles  were  computed  along  the  ridges  indicated  in  Figure  2  and  are  shown  in 
Figure  3. 


Theory 

To  estimate  the  expectation  of  holes  of  various  heights  existing  in  the  ridges,  the  following 
calculation  was  carried  out.  First,  the  random  function  W(x)  defined  by  //(x)  height  of  pressure 
ridge  with  the  variable  x  being  measured  along  the  length  of  the  pressure  ridge  was  introduced.  We 
assumed  that  //(x)  is  normally  distributed  and  homogeneous.  These  assumptions  are  reasonable  in 
view  of  the  variation  of  heights  along  the  ridges  shown  in  Figure  3. 


Figure  3.  Digitized  profiles  at  5-m  intervals  of  ridges  in  Figure  2. 


Figure  3  (coat'd). 
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The  mean  height  and  autocovariance  functions  were  defined  by 


<  //(<)  >  =  H  (2) 

Cfx,,  x.)  =  <  [H(Xj)  -  H][//(xj)  -  H]  >  (3) 

where  brackets  denote  expectation  values  and  C(x.,  x j )  =  C ( j x ^  -  x.|)  since  H  is  homogeneous. 
Following  Smith  and  Nakano  (1973),  if  the  covariance  matrix  ytj  =  C(x,,  x.)  is  nonsingular,  then 
the  multivariate  probability  density  function  ln(H  v  H ...  //n.  x,.  xa  ...  xj)  (probability  of  point 
*j  being  at  a  height  H  with  xa  at  height  H2,  etc.)  is"  given  hy 


'n  - 


|A|‘ 


(2  «)B 


-  exp  -  _  A  .  (H.  -  H)(H.  -  H) 

n/2  p  1J  1  I 


(4) 


where  repeated  indices  are  summed  over 


Aij  "  (Fijrl  (5) 

and  IA|  -  determinant  of  A^. 

Using  the  multivariate  pdf.  the  probability  P  that  there  is  no  point  5  units  above  the  mean 
height  //  in  a  passage  of  width  nAx  where  Ax  is  spacing  between  sample  points  is  given  by 

PlHiXj)  <  H  +  S;  H(Xjj)  <  H  *  5;  ...| 


//+5 
=  f 


lll8 
/  dH 


dH n  «»i 


Hn-*l 


Xn)- 


(6) 


To  carry  out  the  calculation  indicated  in  eq  6,  it  was  necessary  to  determine  values  of  the  auto¬ 
covariance  function  C(Xj,  xp  for  typical  ridge  profiles. 


Results 

Using  the  digitized  profiles  shown  in  Figure  3,  autocovariance  functions  were  calculated  and 
compared.  A  summary  of  the  mean  heights  and  variances  for  ridges  2-9  is  given  in  Table  I.  Ridge 
1  was  not  included  in  the  calculation  because  it  was  of  insufficient  length  for  statistical  purposes. 
In  general,  it  was  found  that  the  autocorrelation  functions  [C(r)  normalized  to  1  at  r  =  01  of  the 
various  ridges  were  quite  similar.  For  an  average  probability  of  finding  ridge  passage,  the  C(r) 
va’ues  for  all  ridges  used  were  averaged  at  each  value  of  r  (Table  II)  and  the  integral  in  eq  6  was 
performed  for  holes  up  to  25  m  in  width.  The  results  are  shown  in  Figure  4.  Since  the  auto¬ 
correlation  factor  rather  than  the  autocovariance  factor  was  used,  all  heights  were  normalized  to 
the  square  roots  of  the  respective  ridge  height  variances.  The  average  variance  was  used  to 
calculate  the  ridge  passage  probabilities.  In  addition  to  the  probabilities  obtained  from  the  average 
autocorrelation  values,  the  limiting  cases,  if  there  were  no  correlation  at  all  between  the  adjacent 
points  on  the  ridge,  are  shown  in  Figure  4.  These  represent  the  worst  possible  cases  of  the 
probability  of  finding  passage  based  on  the  averaged  variance  of  the  ridges. 

To  use  Figure  4  to  calculate  the  probability  of  finding  at  least  one  "hole”  in  a  given  ridge, 
an  encountered  ridge  can  be  broken  into  adjacent  sections  of  length  equal  to  the  desired  hole 


P'OfiOb'lity  b<  Finding  Passage 
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Passoge  Width,  m 


Figure  4.  Average  probability  of  finding 
passage  vs  passage  width.  The  h0  values 
refer  to  height  of  passage  relative  to  the 
wean  ridge  height;  i.e.,  h0  =  -0.67  refers 
to  holes  0.67  m  below  the  wean  ridge 
height.  The  large  dashed  curves  were 
calculated  assuming  no  correlation  be¬ 
tween  adjacent  points  on  the  ridge  and 
the  solid  curve  was  calculated  from  the 
average  autocorrelation  function.  The 
data  points  indicate  the  probabilities  of 
finding  passage  ohtained  by  actually 
sampling  the  ridges.  The  points  are  the 
mean  value  of  the  probability  with  error 
b;vs  given  by  the  standard  error  of  the 
mean  (or  the  seven  ridges  sampled. 


Table  I.  Mean  heights  and  variances 
for  ridges  2-9.* 


Ridge 

ofm2l 

H  (m) 

Length  (m) 

2 

0.43 

2.03 

885 

O 

u 

0.38 

1.83 

930 

4 

0.32 

1.66 

1200 

5 

0.33 

1.69 

926 

6 

0.47 

1.91 

1410 

7 

0.32 

1.90 

1145 

:} 

Avg 

0.87 

0.45 

2.09 

1.87 

1102 

1103 

*  Ridges  8  and  9  are  two  parts  of  the  same 
ridge. 


Table  II.  Average  autocorrelation  function 
C(r)  for  ridges  2-9. 

1  lag  =  5  m. 


Lag  no. 

C(  o 

0 

1.0 

1 

0,508192 

') 

0.348886 

3 

0.294383 

4 

0.234582 

5 

0.192785 

6 

0.187736 

7 

0.187511 

8 

0.185677 
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a.  50-m  segment  p“r-  b.  50-m  segment  per¬ 
centage  with  holes  ot  rent  age  with  holes 

height  1  m  below  for  of  height  equal  to  for 

lower  than)  the  mean  lower  than)  the  mean 

ridge  height.  ridge  height. 

Figure  a.  Percentage  of  50-m  ridge  segments  containing 
holes  of  various  widths.  Total  number  of  50-m  segments 
examined  was  172.  The  percentage  of  segments  with  at 
least  one  hole  up  to  20  m  in  width  was  100 %  for  hole 
heights  1  m  above  the  mean  ridge  height. 

width  which,  to  a  first  approximation,  may  be  considered  independent.  The  probability  of  a  section’s 
being  a  hole  is  given  in  Figure  4;  if  the  vehicle  finds  no  hole  then  the  next  section  may  be  looked 
at.  Figure  4  shows  that  if  the  probability  of  an  acceptable  hole  of  length  l  is  P  then  the  prob¬ 
ability  of  finding  at  least  one  hole  in  a  length,  say  lOt,  is  1  -  (1  -  P)10. 

The  results  in  Figure  4  generally  indicate  that  the  probability  of  finding  wide  holes,  both  from 
the  theoretical  curves  and  from  the  data  points,  is  not  good.  This  is  apparent  from  the  way  the 
curves  fall  off  and  their  good  agreement  with  the  mean  value  obtained  by  sampling  the  profile  as 
indicated  by  the  data  points.  As  another  check  on  this  conclusion,  the  eight  ridge  profiles  were 
broken  into  50-m  segments  and  the  number  of  segments  with  holes  of  different  widths  catalogued. 

The  results,  as  shown  in  Figure  5,  generally  verify  the  conclusion  that  there  are  few  wide  holes  in 
ridges.  Note  that  Figures  4  and  5  are  not  directly  comparable  since  Figure  5  shows  the  number  of 
50-m  segments  having  holes  of  a  given  width,  whereas  Figure  4  gives  the  average  probability  of  a 
given  number  of  random  selected  adjacent  points  being  below  a  given  height. 

Conclusions 

Based  on  the  results  shown  in  Figures  4  and  5,  vehicle  trafficability  can  broadly  be  classified 
in  relationship  to  the  terrain  in  the  following  manner.  Vehicles,  in  their  ability  to  negotiate  ridges. 
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fall  into  two  classes.  If  a  vehicle  is  able  to  negotiate  only  small  heights  (below  the  mean  ridge 
eight;,  it  must  be  quite  narrow  since  the  availability  of  holes  drops  off  rapidly  with  the  width 
necessary  for  passage.  This  means  that  a  vehicle  like  a  snowmobile  with  good  maneuverability 
could  probably  negotiate  the  pack  ice  although  considerable  time  for  hole-searching  would  be 
necessary.  If  a  vehicle  has  considerable  beam  width  (  10-15  111).  the  probability  of  finding  a  hole 
is  so  low  that  the  vehicle  must  go  mound  the  ridges.  The  ramifications  of  this  vis-A-vis  .1  more 
complete  trafficability  analysis  are  discussed  by  llibler  and  Ackley  (1973).  The  general  nature 
0  t  ie  complete  trafficability  analysis  is  what  one  would  intuitively  expect:  since  ridges  often 
extend  beyond  2.5  km,  the  presence  of  on-t lie-average  more  than  0.5  impassable  ridges/km  makes 
it  difficult  for  the  vehicle  to  proceed  rapidly  and  economically.  This  may  be  translated  into  the 
necessary  height  clearance  ability  of  a  vehicle  by  using  ridge  statistics  (llibler  et  al.  1972) 
compiled  for  different  regions  of  the  Arctic  during  different  seasons  (Mock  et  al.  1972). 
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